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ABSTRACT 

The pressure disturbances generated by an instabil- 
ity wave in the shear layer of a supersonic jet are 
studied for an axisymmetrie jet inside a lined cylin- 
drical duct. For the supersonic jet, locally linear sta- 
bility analysis with duct wall boundary conditions is 
used to calculate the eigenvalues and the eigenfunc- 
tions. These values are used to determine the growt h 
rates and phase velocities of the instability waves 
and the radial pressure disturbance patterns. The 
study is confined to the dominant Kelvin-Helmholtz 
instability mode and to the region just downstream 
of the nozzle exit where the shear layer is growing 
but is still small in size compared to the radius of 
the duct. Numerical results are used to study the 
effects of changes in the outer How, growth in the 
shear layer thickness, wall distance, wall impedance, 
and frequency. Results indicate that the 1 effects of 
the duct wall on shear layer growth rates diminish 
as the outer How increases. Also, wall reflections 
cause variations in growth rates depending on wall 
height and Strouhal number. These variations are 
due to the phase relationship between the outgoing 
and the reflected incoming pressure disturbances at 
the shear layer. The growth rate variations can be 
reduced and the maximum growth rate minimized 
by keeping the imaginary part of the impedance neg- 
ative'. 

INTRODUCTION 

A confining wall affects the supersonic jet mix- 
ing layers that arc 1 associated with combustion pro- 
cesses, 1 engine tests conducted inside* an enclosed 
facility, 2 and jet noise reduction efforts/* In the cir- 
cular configuration considered for this paper, the su- 
personic jet issues from a nozzle* that has its axis co- 
incident with the centerline of a surrounding cylin- 
drical duct. Between the jet and the* confining wall 
of the duct, a lower speed outer How exists. The* 
shear layer between the two flows is initially thin 
and unstable. This instability grows like a wave as 
the How moves downstream. As a result, this insta- 
bility wave generates perturbations that propagate* 


away from the shear layer. Those perturbations that 
propagate toward the wall can reflect off it and re- 
turn toward the shear layer to interact and modify 
the growing instability in the shear layer. Thus the 
presence of the wall, both if the wall is hard and if it 
lias an acoustic lining, affects the characteristics of 
the instability waves that govern both the mixing in 
the high speed shear layer and the noise generation. 

Stability analysis was used by Tam and Hu' 1 to 
study the supersonic mixing layer inside a rect- 
angular channel with hard walls. In addition to 
the Kelvin-Helmholtz (KH) instabilities, they found 
from using a vortex sheet model that the coupling of 
the unsteady motion with acoustic modes reflected 
from the walls resulted in a new family of supersonic 
instabilities with supersonic convective Mach num- 
ber. The analysis also included the effects of shear 
layer thickness on the spatial growth rates of t wo- 
und three-dimensional supersonic instability waves. 
Zhuang et al/’ showed for a constant shear layer 
thickness that the walls most affected the shear layer 
instability when the convective Mach number was 
greater than one. As the distance between the walls 
decreased, the maximum amplification rate reached 
a maximum value and then decreased. The stability 
characteristics for subsonic convective Mach number 
waves were unaffected by wall height. Others have 
found similar results for shear layers in a rectangular 
channel . 1(< 

The stability analysis has been performed for a 
cylindrical supersonic vortex sheet inside a cylindri- 
cal duct with hard walls. 2,7-9 ViswanathaiT showed 
how the KH instabilities for a Mj = 4.0 jet with no 
outer flow were affected by the wall height. The ef- 
fect appeared as oscillations in the growth rates due 
to interaction between the internal Mach wave sys- 
tem of the jet and the reflected waves from the wall 
as the wall moved closer to the vortex sheet,. Fur- 
thermore, the KH mode had growth rates typically 
higher than the supersonic instabilities. Chang and 
Kuo 8 found similar results for a Mj = 4.5 jet with 
two different outer Hows. M a ~ 4.0G and M n = 1.56. 
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It was also shown that the axisymmetric I\H mode 
growth rates decreased as the outer flow velocity in- 
creased relative to the inner flow velocity. When 
shear layer thickness was included in the analysis, 
the growth rates for the instabilities were diminished 
and the presence of the hard wall had a greater af- 
fect on I\H mode growth rates than on supersonic 
mode growth rates. 7 

Hu : * recently studied the effects of two paral- 
lel walls with sound absorbing lining on a two- 
dimensional jet confined inside a rectangular duct. 
The analysis was performed for two-dimensional in- 
stability waves with mean flow conditions Mj = 2.0, 
M 0 — 0.2, and Tj/T ( , = 4.0. Using a vortex sheet 
model, the lined walls were shown to be most ef- 
fective when the phase velocity was supersonic rela- 
tive to the outer flow and when the pressure eigen- 
function showed large perturbations near the wall. 
In general, growth rates were reduced when the 
wall impedance was finite. A continuous mean flow 
model based on a hyperbolic tangent function was 
used to include the effects of shear layer thickness. 
For the instability wave with the largest growth 
rates, the increasing shear layer thickness resulted 
in lower growth rates. Also, the presence of a lined 
wall was more effective in reducing growth rates than 
a hard wall. 

In this paper, we extend the study of the effects of 
a confining wall with sound absorbing lining on shear 
layer instability waves to the cylindrical duct case. 
Many variables are available to study including the 
inner jet stream velocity and temperature, the outer 
stream velocity and temperature, the axial location, 
the jet shear layer width, the duct wall height, the 
wall impedance both real and imaginary parts, and 
the frequency. By necessity, we must limit our pre- 
sentation of results. The jet Mach number is fixed 
at 2.0 with a temperature ratio Tj/T () — 2.0. We 
are primarily interested in the application of acous- 
tic wall treatment in an ejector designed to reduce 
supersonic jet noise. Inside the ejector, a shear layer 
between the two flow streams exists only for a short 
distance downstream of the nozzle exit before the 
tw'o flows more fully mix together. Consequently, 
the results presented herein will be confined to small 
shear layer widths near the nozzle. Thus, this study 
concentrates on the effects of w r all treatment on the 
initial growth characteristics of the shear layer. We 
will only consider the Kelvin-Helmholtz type mode 
in this study since it is the dominant instability for 
both the free and confined jets at the above flow 
conditions. The next section describes the formu- 
lation of the instability wave model using a finite 
impedance boundary condition. The model is used 
to determine the growth rates and phase velocities 
of the instability waves that govern the pressure dis- 
turbance pattern generated by the jet shear layer. 
This is followed by a presentation and a discussion 
of the numerical results. 


INSTABILITY WAVE FORMULATION 


A cross-section of the model configuration for a 
supersonic axisymmetric jet inside a. cylindrical duct 
is shown in Figure 1. The high speed flow emerging 
from the nozzle with velocity Uj and static temper- 
ature Tj has an initially thin shear layer that is in- 
herently unstable even in the absence of viscosity. 
An instability wave begins to grow rapidly and con- 
tinues to grow* at reduced rates as the shear layer 
spreads. During this process, pressure disturbances 
created by the instability propagate away from the 
shear layer toward the wall, reflect off the wall, and 
return to interact with and modify the instability 
wave in the shear layer. This process is assumed to 
be governed by the linearized, in viscid, compressible 
equations of motion. 


Disturbance Equations 


For a supersonic jet, the shear layer grow r s slowly 
in the axial direction. This allows a locally parallel 
flow 7 approximation to be used to obtain the solution 
for the disturbance quantities. 4 The fluctuating dis- 
turbances are represented as waves traveling through 
a nonuniform medium. As an example, the pressure 
disturbances are given by 


;/(r, 0. j\ t) — p(r 1 x) 
x exp \i 



( 1 ) 


where p(r, x) represents the radial distribution of the 
pressure disturbance at each axial location, a(x) is 
the local complex w r ave number (a = a r -b m, and 
—a, is the local growth rate), 0 is the azimuthal an- 
gle, n is the mode number, and exp(— iuit) is the 
harmonic time dependence. The other disturbance 
quantities have 4 the same form. The linearized equa- 
tions governing the disturbances can be combined 
to obtain a single equation called the compressible 




u r T i' M i 
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Figure 1: Schematic diagram in the x-r plane of ail 
axisymmetric jet flowing into a cylindrical duct with 
coflowing outer stream. 
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P = 0. (2) 


This equation lias been nondimensionalized as fol- 
lows: spatial coordinates by i? ; . time and radial fre- 
quency by Uj/Rj , velocity by l/j, density by pj, and 
pressure by pj Uj where the subscript j indicates jet 
exit conditions. To solve this equation, the radial 
variation of the mean flow velocity u and density p 
must be known quantities provided from either an 
analytic expression or a numerical calculation. 

The general solution to equation (2) at any axial 
location is written 


p(r) = ACf (r) + B(*(r) (3) 

using unknown functions with radial variations. 
Outside the jet, the mean flow becomes uniform and 
equation (2) takes the form of Bessel's equation. The 
solution can be readily written as 

p„ = A 0 Hl l) (iX(a)r) + B 0 H™ (i\(a)r) (4) 

where 

A(n) = [o- - - rt1I„) 2 ] ' /_ . (5) 


H { n l) is a Hankel function of t.he first kind. It rep- 
resents a wave outgoing from the shear layer toward 
the wall. H\i ] is a Hankel function of the second 
kind and represents the wave reflected from the wall. 
As r — > 0, equation (3) must be finite. 

At the wall, the continuity of particle displace- 
ment is used as the boundary condition.** The kine- 
matic equations are combined in terms of pressure 
to get 


- au 0 )~p 0 + — 


iuj Z* dp () 


PoPjUj 0r 


= 0 


( 6 ) 


where Z* is the dimensional acoustic impedance. 
The acoustic impedance is usually referenced to a 
typical density and speed of sound. We use ambient 
conditions for density p a and speed of sound c n to 


get 


^ _ y /=— 1 

rr ” 


f r — v r a » r ( ^ ) 

J j 

where Z = Z* / p a c„ is the specific acoustic 
impedance of the wall treatment. This formulation 
for the wall treatment condition is based on assum- 
ing a locally reacting boundary. Substituting equa- 
tion (4) into equation (G) and evaluating the result- 
ing equation at the radial location h of the wall for 
the unknown B a < we get the equation 


Po = Ao [Hi! HiXr) - C(h)H™(iXr)\ (8) 


C(h) = 


(^-<m n y 2 Hl"(rXh) 
(*> - (iXh) 


;X f^ Hl, u '(iXh) 

P< > 


~TJT~ 


H^'(iXh) 


( 9 ) 

is obtained after substituting B (l back into equation 
(4). As // — > oc. equation (8) reduces to the free jet 
case 10 and as Z — » oc, the wall becomes acoustically 
hard and equation (9) reduces to the appropriate 

form (i\h) / H { u 2)f {i\h) as found in Ref. [9]. 

Equation (2) and its boundary conditions create 
an eigenvalue 1 problem with eigenvalue a. The cen- 
terline boundary condition at r — 0 depends on 
whether the pressure eigenfunction is either axisym- 
met.ric, n — 0, or non-axisymmetric, n ^ 0. In oqua- 
tion form, this is written as 


& 

P 


= 0, n = 0 
= 0 , v ± 0 . 


(10) 


A finite difference approximation is used to discretize 
equation (2). To incorporate 1 the outer boundary 
condition, the constant A a must be eliminated from 
equation (8). This is done by using a ratio of the 
disturbance pressures at the two outer grid points, 

PK + 1 W.VW.v+i) ~ C(h)H:r ] (,\r x+l ) 

Pft _ Hi, l) {iXr N )-C[h)Hl?\iXr N ) 

where N is the number of grid spacings. The 1 lo- 
cal eigenvalue is found from the resulting diagonal 
matrix using a Newton-Raphson iteration for refine- 
ment. The local eigenfunction, p(r,.r). is then de- 
termined within a constant using the inverse power 
method. During this latter calculation, the eigen- 
function is normalized at the point where the axial 
mean flow velocity gradient, is a maximum in the 
radial direction. 


Mean Flow 

In order to complete the stability analysis, the 
mean How properties are needed. Often equations 
consisting of simple analytic functions are used to 
describe the mean flow. However, in order to eas- 
ily account for the effects of velocity and tempera- 
ture differences across the supersonic jet shear layer 
on shear layer spreading, a numerically generated 
mean flow is used. The jet is assumed to be per- 
fectly expanded and for the region of interest near 
the nozzle, the static pressure is assumed to be con- 
stant. The set of compressible, Reynolds averaged, 
boundary layer equations with a modified turbulence 
model is used as the basis for the numerical mean 
flow analysis. The details of the equation develop- 
ment. the numerical solution, and comparisons of 
calculated results to measured mean flow data are 
given in Ref. [11]. 
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Figure 2: The integrated instability wave amplitudes 
to the free jet shear layer half-width b = 0.2. Modes: 

( ) n = 0; (- - - -) n = 1; ( ) n = 2. 

Outer flow: thick lines, AI a = 0.03; thin lines, M a = 
0.6; peaks shewn for M 0 = 0.1 to 0.5. 


NUMERICAL RESULTS 

The' calculations were performed for a jet with exit 
Mach number Mj — 2.0. The jet was hot relative 
to the outer stream with Tj/T 0 — 2.0. The outer 
stream had the same static temperature as the am- 
bient conditions for the source of the outer stream. 
The velocity ratio between the outer stream and the 
inner jet stream was varied from 0.01 to 0.21 result- 
ing in outer stream Mach numbers M 0 of 0.03 and 
0.1 to 0.6 in increments of 0.1 for a total of seven 
operating conditions. 

In addition to the flow conditions, the frequency u 
in equation (2) must be given. A study of the insta- 
bility wave amplitudes for a free jet shear layer cal- 
culated from the nozzle to the downstream location 
where the shear layer half-width b is 0.2 shows peak 
amplitudes for the n = 1 and n = 2 helical modes 
in the Strouhai number range (St — u;/7r = JDj/Uj 
where Dj = 2 Rj) between 0.2 and 0.3, see Figure 2. 
As the outer flow Mach number increases, the peak 
wave amplitudes increase and shift slightly to higher 
Strouhai numbers. Thus, the calculated results will 
be shown for Strouhai numbers in this peak region. 
In addition, the results presented in this paper will 
be for the n = 1 and n = 2 instability wave modes 
since they are clearly more dominant than the n — 0 
axisv inmetric instability wave mode. 

Hard Wall 

The instability wave growth rates near the noz- 
zle exit are first examined for the hard wall case 
as a function of the Strouhai number St , the duct 
wall height /i, and the outer stream Mach number 
M (> . Figure 3 shows growth rate contours at several 



(c) (d) 



Figure 3: Local growth rate contours as a function 
of h and St for the n = 1 mode of an Mj = 2.0 jet in 
a hard wall cylindrical duct, b = 0.05. Outer flow: 
M 0 = (a) 0.03; (b) 0.1; (c) 0.3; (d) 0.6. Phase veloc- 
ities supersonic relative to outer flow and subsonic 
relative to inner flow except inside thick dashed line 
(subsonic re. outer flow) and inside thin dashed line 
(supersonic re. inner flow). 


outer stream Mach numbers for the n = 1 mode as a 
function of wall height and Strouhai number at the 
axial location where the shear layer half-width b is 
0.05. For M a = 0.03, the growth rates oscillate as 
both the wall height and Strouhai number change 
becoming large near the lower wall height and the 
largest Strouhai number. This behavior is similar 
to those growth rate characteristics calculated for a 
cylindrical vortex sheet. 7 These growth rate oscilla- 
tions are due to the location of the wall relative to 
the standing wave like pattern resulting from reflect- 
ing pressure disturbances interfering with the pres- 
sure disturbances generated by the shear layer insta- 
bility. As an example, Figure 4(a) shows the changes 
in the radial distributions of the wave modes for the 
local maximum and minimum growth rates as the 
wall height changes for St — 0.45 in Figure 3(a). To 
indicate the relative amplitude of these modes, the 
radial distributions are multiplied by their respec- 
tive local growth rates and then normalized by the 
growth rate at h — 1.65. Local nodes and antin- 
odes occur that are due to the interference between 
outgoing and incoming reflected disturbances. The 
total pressure disturbance amplitude is always a lo- 
cal maximum at a hard wall. As the reflected dis- 
turbance interacts with the shear layer, the growth 
rate can be enhanced or reduced. Figure 4(b) shows 
the outgoing and the incoming pressure disturbance 
patterns separately (note that the pressure distur- 
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bailee reflects in phase at the wall). The major dif- 
ference between the instability wave having a local 
maximum growth rate or a local minimum growth 
rate is the phase relationship between the outgoing 
and incoming disturbances. For a local maximum 
growth rate, the phase differences between the in- 
coming and the outgoing disturbances just outside 
the shear layer are -148.5°, -152.6°, and -152.6° 
when the wall heights are 1.65. 2.70, and 3.75, re- 



Figure 4: Radial distribution of wave modes for 
conditions in Figure 3(a), St = 0.45. (a) Pres- 

sure magnitudes for different, wall heights normal- 
ized bv «,•(/*. = 1.65) = -1.1762. ( ) h = 

1.65: ( ) h = 2.33; ( ) h = 2.70; 

( ) h = 3.38: ( ) h = 3.75. (b) Real 

part of pressure disturbance, (■"^■"■) outgoing; 
( ) incoming. 



Figure 5: Radial distribution of n = 1 wave modes 
at maximum growth rate with St — 0.45 for an 
Mj = 2.0 jet in a hard wall cylindrical duct, h — 
0.05. Normalized by (\j(M a = 0.03) = —1.1762. 

( ) M 0 = 0.03, h = 1.65; ( ) M a = 

0.1, h = 1.65: ( ) M 0 = 0.2, h = 1.73: 

( ) M 0 = 0.3. h = 1.80; (----) M a = 

0.4. h = 1.88; (— mm mm) M a = 0.5, h = 2.03; 
(mm - — ) M a = 0.6. h = 2.10. 


spectively. The ])hase differences for a local mini- 
mum growth rate are 52.4° for h ~ 2.33 and 66.7° 
for h = 3.38. Thus, the pressure disturbances are 
more out of phase for a local maximum growth rate 
than for a local minimum growth rate. 

As M a increases, the variations in the growth rate 
contours in Figures 3(b) to (d) decrease and the gen- 
eral trend is toward lower maximum growth rates 
for St < 0.45. This latter trend was observed in 
calculated results for both free jets with a coflow- 
ing stream 12 and confined cylindrical vortex sheets. 8 
The largest growth rate within the calculated region 
for each of the four flow conditions in Figure 3 oc- 
curs at St ~ 0.45 and at a lower wall height. The 
radial distributions for these growth rates are shown 
in Figure 5 along with the largest, growth rate results 
for the other three flow conditions. The relative de- 
crease in the amplitude of the mode as M v increases 
is shown by normalizing the radial mode distribution 
using the local growth rate. The radial distribution 
is multiplied by its local growth rate and then nor- 
malized by the local growth rate at M 0 — 0.03. The 
decrease in amplitude is due to the change in the 
direction of the disturbance wave vector as the flow 
increases. From equation (1), we (‘an define a phase 
function for a disturbance wave 

<f> = S(r) -»'t (12) 


where 

.r 

S(r) = if n (iXr) + I o r (\)r/\ + nO (13) 
b 

and ip n (i\r) is the real phase function for H ( n l) (i\r) 
in the case of an outgoing disturbance wave. To fol- 
low the local wave front, $ in equation (12) must 
remain constant and following Goldstein. l '* an in- 
cremental change in <t> to lowest order is 

V S • dr = Ljdt. (14) 


From this equation, the unit normal to the wave 
front is defined by VS/|VS| with components 


os = o^ 

[ J '- dr dr 

[VS1 =I^ = !i 
e ^9 r go 

r^ci 0S 

[vsl = aF = n - 


and with direction cosines 


COS0 r — 

cos (j)# = 

COS0j. = 


1 d<p„ 
|VS| dr 
1 n 


|V5| r 
o, 

|VS| 


(13) 


( 16 ) 
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Using these results, we determine from the eigen- 
function results in Figure 5 that the wave front di- 
rection <j) x changes on average from 57.00° ± 0.05° 
to 42.75° db 0.38° as the outer flow Mach number 
increases from 0.03 to 0.6. The lower angle means 
the wave vector points more towards the :r-axis in 
the downstream direction and less of the reflecting 
disturbance returns to the shear layer to have any 
effect. As a consequence, the growth rate and the 
mode shape are approaching that of the free jet as 
M 0 increases. 

Finally, the regions of supersonic or subsonic rel- 
ative phase velocity are demarcated by dashed lines 
in the plots of Figure 3. For a phase velocity super- 
sonic relative to the inner stream jet velocity, 

1 1 ~ ('ph | > 1 00 

and for a phase velocity supersonic relative to the 
outer stream, 

|«o -t-ph I > -rr-. (18) 

Mo 

where Cph — uj/a v is the local phase velocity of the 
instability wave. In Figures 3(a) and (b), the phase 
velocities are all supersonic relative to the outer flow. 
A small portion of the growth rate plot becomes sub- 
sonic relative to the outer flow at M 0 = 0.3 (lower 
left of part (c)). This region increases in size as 
M 0 increases as seen in part (d). In the small re- 
gions outlined by the thin dashed lines in parts (a) 
to (c), the phase velocities are supersonic relative 
to the inner flow. Most of the growth rates lie in 
the region where the phase velocity is subsonic rel- 
ative to the inner flow. The most important results 
for the wall to affect the shear layer growth rates 
is that the phase velocity be supersonic relative to 
the outer flow, especially when the wall has acoustic 
lining/ The flow conditions chosen meet the crite- 
rion of having the phase velocity supersonic relative* 
to the outer flow when lined walls are present. 

Lined Wall 

It has been previously shown that instability wave 
growth rates decreased as the real part, of the wall 
impedance decreased toward one for a fixed imagi- 
nary part and wall height/ We confirmed this be- 
havior for the cylindrical duct case with similar op- 
erating conditions and impedance boundary condi- 
tions. In this paper, we present results for a range 
of wall impedances that are found in many practi- 
cal applications. The real part of the impedance is 
always positive with values chosen to be between 0 
and 2. The imaginary part is either positive or nega- 
tive with values between —2 and 2 chosen for study. 
These impedance* values are realizable with practical 
sound absorbing materials. 


Growth Rates versus Impedance 

The initial growth rates of the n — 1 shear layer 
instability wave modes were calculated for a wall 
height h — 2 and a Strouhal number St — 0.2. With 
seven outer flow conditions, a large number of calcu- 
lations were performed over the designated region in 
the impedance plane. To compare this large number 
of calculated results, we chose to compare* growth 
rates affected by the presence of a wall impedance to 
the growth rate of a free jet shear layer. The growth 
rates were mapped on to the impedance plane and 
contours were drawn through the data correspond- 

(a) (b) 





Figure 6: Local growth rates for the n — 1 mode 
of an Mj — 2.0 jet in a lined cylindrical duct, 
h . = 2, St = 0.2. Contours in the impedance plane 
for: (a) b « 0.0; (b) b = 0.05; (c) b = 0.10; 
(d) b = 0.15. Contours for free jet growth rate 

at M () = 0.03 ( ) and M a — 0.1 to 0.6 

( ) following arrow. Arrow points from re- 

gion where growth rates lower than free jet to region 
where growth rates higher than free jet. 
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ing to the 4 free jet growth rate; thus, marking the 
boundary between regions where 4 the growth rate 
was higher than the free jet growth rate and regions 
where it was lower. Examples are shown in Figure G. 
Eac h part of the figure 4 shows seven contours corre- 
sponding to the free jet growth rates for the 4 seven 
outer flow conditions at a fixed axial location deter- 
mined by the jet shear layer half-width. The arrows 
point in the direction that the* contour moves with 
increasing outer flow Mach number M n and point 
into the region of higher growth rates. As M a in- 
creases, the region of higher growth rates decreases 
in size and as the shear layer half-width increases, 
this area expands especially at lower M (> . 

Figure G only shows the 4 separate regions in the 
impedance plane 4 between regions with growth rates 
higher or lower than that of the free 4 jet. Figure 7 
shows the range of those growth rates. For each 
shear layer half-width location, the figure 4 shows the 
largest growth rate, the smallest growth rate 4 , and 
the free jet growth rate within that calculated re- 
gion of the impedance plane. As b increases, the 
growth rates generally decrease 4 and as M a increases, 
the growth rates decrease and the range e)f increased 
growth rates, the difference between the largest, 
growth rate and the free 4 jet growth rate 4 , also de- 
creases. The largest growth rates shown in Figure 7 
lie 4 along the imaginary axis e>f the impedance 4 planes 
in Figure G ne 4 ar the center of a curie that approxi- 
mate's the shape 4 of the free jet grow th rate 4 contours. 

Figure 8 shows how the fire jet growth rate 4 con- 
tours change in the impedance plane wdien the 4 duct 
wall height is changed with b = 0.05. St — 0.2, and 
n = 1. The 4 duct wall height h starts at LG and 
progresses to 2.0, 2.5, and 3.0. At h = 1.6, the* 
contours change location a small amount. As h in- 
creases, the 4 contours unfold creating larger regions 
of higher growth rates. The contours at lower h are 4 
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M o 


Figure 7: Range of local growth rates for condi- 
tions in Figure G. ( ) free jet growth rate; 

(■■ — mm) largest growth rate 4 ; smallest 

growth rate, (a) b ~ 0.0; (b) b — 0.05, thin lines 
correspond to conditions in Figure 8; (c) b — 0.01; 
(d) b = 0.15. 


circular in shape, wrapping around the 4 point of the 4 
largest growth rate 4 . The 4 contours unfold as h in- 
creases until they start to w r rap around the point of 
the smallest growth rate. This alse> may be viewed as 
a movement of the largest growth rate 4 point toward 
larger positive imaginary impedance values, outside 
the impedance range 4 shown, as h increases and a 
movement of the smallest growth rate 4 point into the 
region from more negative? imaginary impedance 4 val- 
ues. There may be more local maximum and mini- 
mum growth rate 4 points that lie outside this region 
of the impedance plane 4 . Even though the region of 
increases! growth rates get larger as h increases, the 


(a) (b) 




Figure 8; Local growth rates for the n — 1 mode 
of ail Mj = 2.0 jet in a lined cylindrical duct, 
b - 0.05. St — 0.2. Contours in the impedance 4 
plane for: (a) h — 1.6; (b) h = 2.0; (c) h = 2.5; 
(d) h = 3.0. Contours for free jet growth rate 

at M 0 = 0.03 ( ) and M t) = 0.1 to 0.G 

( ) following arrow. Arrow points from re- 

gion where 4 growth rates lenver than free jet to region 
wtiere 4 growth rate 4 s higher than free 4 jet.. 
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range of those growth rates decrease as shown in 
Figure 7(b). The increasing outer flow decreases the 
range even more until for the M 0 = 0.6 growth rates 
in the impedance plane of Figure 8(d), there is only 
a 10% change from smallest to largest growth rates 
at h = 3.0. 

The results for the n = 2 helical mode at St = 0.2 
follow the same trends as those seen for the ?i = 1 
mode above though with generally higher growth 
rates. An example of the calculated growth rate re- 
sults in the impedance plane is shown in Figure 9 
in terms again of the free jet growth rate contours. 
These plots for b = 0.05 and h — 2 should be com- 
pared to the n = 1 results shown in Figures 6(b) 
and 7(b). 

Calculations were also performed at a Strouhal 
number of 0.3, a duct wall height of 2, and a shear- 
layer half-width of 0.05. The free jet growth rate 
contours are shown in Figures 10(a) and (b) for the 
n = 1 and n = 2 modes, respectively. The corre- 
sponding growth rate ranges as a function of M a are 
shown in parts (c) and (d). Note the similarity be- 
tween Figures 10(a) and 8(c). The duct wall heights 
and Strouhal numbers for these two cases are such 
that St(h — 1) is a constant. This indicates that an 
important length scale near the nozzle may be the jet 
radius to duct wall distance. However, the 1 growth 
raters may be higher and the range of growth rates 
is certainly larger for the St = 0.3 case than for the 
St — 0.2 case, compare Figure 10(c) to Figure 7(b) 
with h = 2.5. 


(a) (b) 



m q m o 


Figure 10: Local growth rates of an Mj = 2.0 
jet in a lined cylindrical duct, h — 2. Si — 0.3. 
b = 0.05. Contours in the impedance plane: (a) 
n = 1; (b) n — 2. Contours for free jet growth 

rate at M 0 = 0.03 ( ) and M 0 = 0.1 to 0.6 

( ) following arrow. Arrow points from re- 

gion where growth rates lower than free jet to re- 
gion where growth rates higher than free jet. Range 
of local growth rates: (c) n = 1; (d) n = 2. 

free jet growth rate; — — i largest 
growth rate; (««■—) smallest growth rate. 



Figure 9: Local growth rates for the n = 2 mode 
of an Mj = 2.0 jet in a lined cylindrical duct, 
h — 2, St = 0.2, b = 0.05. (a) Contours in the 
impedance plane for free jet growth rate at M 0 = 
0.03 ( ) and M 0 = 0.1 to 0.6 ( ) fol- 

lowing arrow. Arrow points from region where 
growth rates lower than free jet to region where 
growth rates higher than free jet. (b) Range of lo- 
cal growth rates. (■—■—) free jet growth rate; 
(— mm —) largest growth rate: (— • - -) smallest 
growth rate. 


Growth Rates versus h and St 


The previous results fixed the geometry and flow 
parameters and explored shear layer growth rates as 
a function of the wall impedance. Next, we select a 
wall impedance and calculate the growth rates as a 
function of wall height and Strouhal number. Using 
the results for M a = 0.3 as an example, Figure 11 
shows growth rate contours for a set of impedance 
values. In the horizontal direction from left to right, 
the real part is 0.4, 1.0, and 1.5 and in the vertical 
direction from bottom to top, the imaginary part 
is —0.8, 0.0, and 0.8. At h — 4.5, the growth rate 
contours on all the plots in the figure are nearly the 
same and mostly vertical indicating that the growth 
rates are about the same as the free jet growth rates 
since the effects of reflections are negligible. As h 
decreases, the effects of the lined wall on growth 
rates first appear at higher Strouhal numbers and 
then progress to lower Strouhal numbers. The ef- 
fects are larger, meaning higher maximum growth 
rates and larger variations in growth rates as the 
wall height changes, when the imaginary part of the 
impedance is positive and for smaller real parts of 
the impedance. 

When the outer flow changes, the previous results 
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Strouhal Number 
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Straubs] Number 
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Figure 11: Local growth rate contours as a function of h and St for the n — 1 mode of an Mj — 2.0 jet in 
a lined cylindrical duct. M 0 = 0.3, b = 0.05. (a) Z = 0.4 4- *0.8; (1>) Z = 1.0 + *0.8; (c) Z = 1.5 + *0.8; (cl) 
Z - 0.4; (e) Z = 1.0; (f ) Z = 1.5; (g) Z = 0.4 - *0.8: (h) Z - 1.0 - *0.8; (i) Z = 1.5 - *0.8. 


have shown that the growth rate changes for each 
wall impedance. If we now fix the wall impedance 
at the extremes of effect in Figure 11, that is Z = 
0.4+ *0.8 where there were large 1 variations in growth 
rates and Z = 1.5 - <0.8 where there were small 
variations in growth rates, the effects of external 
flow on growth rates can be shown as a function 
of wall height and Strouhal number. These results 
are shown in Figures 12 and 13. The four flow con- 
ditions shown in each figure are M (> = 0.03, 0.1, 0.3, 
and 0.6. The M 0 = 0.3 results are the same as shown 
in Figure 11. Thus, we see with lower M (i that there 
are much larger variations in growth rate* and less 


variations as M a increases above 0.3. Tin 1 general 
trends in growth rates discussed for M a = 0.3 in 
Figure 11 also occur at the higher and lower outer 
flow Mach numbers. Figures 12 and 13 can also be 
compared to the hard wall results shown in Figure 3. 
Figure 12 shows generally increased growth rates as 
the w r all height decreases wdien compared to the hard 
wall cases. We find the opposite or generally de- 
creased grow th rates w hen Figure 13 is compared to 
Figure 3. 

CONCLUDING REMARKS 

For a hard w-all confining duct, this analysis con- 
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Figure 12: Local growth rate contours as a. function 
of h and St for the n = 1 mode of an Mj = 2.0 
jet in a lined cylindrical duct with Z = 0.4 4* *0.8, 
b — 0.05. Outer flow: M a = (a) 0.03; (b) 0.1; (c) 
0.3; (d) 0.6. Phase velocities supersonic relative to 
outer flow and subsonic relative to inner flow except 
inside thick dashed line (subsonic re. outer flow) and 
inside thin dashed line (supersonic re. inner flow). 


(a) (b) 




0.06 0.16 0.25 0 36 0.45 0.06 0.15 0.25 O^mT ^46 

Strouhal Number Btrouhal Numb«r 


Figure 13: Local growth rate contours as a function 
of h and St for the n = 1 mode of an Mj = 2.0 
jet in a lined cylindrical duct with Z — 1.5 — *0.8, 
b — 0.05. Outer flow: M 0 = (a) 0.03; (b) 0.1; (c) 
0.3; (d) 0.6. Phase velocities supersonic relative to 
outer flow and subsonic relative to inner flow. 


Armed that the growth rates of Kelvin-Helmholtz 
type instability waves were affected by wall height 
when the waves had phase velocities supersonic rel- 
ative to the outer flow. The results presented here' 
agreed with previous calculated results for the cylin- 
drical vortex sheet case in a hard wall duct where 
the growth rates generally decrease with an increase 
in the outer flow. 8 However, as the wall height and 
Strouhal number change, the growth rate oscillated 
between local minimums and local maximums. The 
range between minimum and maximum decreased 
as the outer flow increased. The phase relationship 
between the outgoing disturbance and the reflected 
incoming disturbance was a factor in determining 
if the growth rate was a local minimum or a local 
maximum. The phase difference was approximately 
— 150° just outside the shear layer for a local maxi- 
mum growth rate and approximately 60° for a local 
minimum growth rate. As the outer flow increased, 
it was shown that the pressure disturbance wave vec- 
tor shifted toward the downstream direction, so that 
for a fixed wall height less reflected disturbance af- 
fected the growth rate. Further analysis using the 
wave vector would reveal direction and level of dis- 
turbance energy. 

The calculated results for the initial growth rates 
for a jet shear layer in a lined duct were presented for 
changes in wall impedance, wall height, outer flow, 
frequency or Strouhal number, and mode number. 
With all conditions fixed except for variations in the 
wall impedance, the growth rate could be made to 
increase above or decrease below the growth rate for 
the free jet. Within the region of the impedance 
plane used in the calculations, a maximum and a 
minimum growth rate existed for the fixed condition. 
This range of growth rates was decreased when the 
outer flow increased at constant shear layer width. 
At fixed flow conditions, an increase in shear layer 
width lowered the growth rates but the range of 
growth rates changed little. 

For a fixed impedance, flow condition, and mode 
number, growth rate variations were calculated as 
the wall height and Strouhal number were varied. In 
general for all the flow conditions studied here, these 
growth rate variations decreased when the imaginary 
part of the impedance was less than zero and in- 
creased when the imaginary part of the impedance 
was greater than zero. With the varied combina- 
tions of wall height, Strouhal number, flow condi- 
tions, and wall impedance, a particular growth rate 
can be higher or lower than the growth rate for a free 
jet or for a jet in a hard wall duct. If it is desirable to 
lower growth rates in general, then those modes and 
frequencies that may tend to grow the most should 
be reduced. Wall impedance values with negative 
imaginary values can be chosen to reduce the initial 
growth rates of these modes. 


NASA TM- 107441 


10 


REFERENCES 


1. D. S. Shin and J. H. Ferziger. Linear stability of 
the confined compressible reacting mixing layer. 
AIAA 31:571 577. 1993. 

2. K. K. Aliuja, K. C. Massey. A. J. Fleming. 
C. K. W. Tam. and R. R. Jones III. Acoustic in- 
teractions between an altitude test facility and 
jet engine plumes theory and experiments. 
Technical Report TR-91-20, AEDC. 1992. 

3. F. Q. Hu. The acoustic and instability waves of 
jets confined inside an acoustically lined rectan- 
gular duct. J. Sound Vib ., 183:841 856, 1995. 

4. C. K. \Y. Tam and F. Q. Hu. The instability 
and acoustic wave modes of supersonic mixing 
layers inside a rectangular channel. J. Fluid 
Meek ., 203:51 76, 1989. 

5. M. Zhuang, P. E. Dimotakis, and T. Kubota. 
The effect of walls on a spatially growing su- 
personic shear layer. Phys. Fluids A , 2:599 604. 

1990. 

6. P. J. Morris and M. G. Giridharan. Tin 1 effect- 
of walls on instability waves in supersonic shear 
layers. Phys. Fluids A, 3:356 358, 1991. 

7. K. Vis wanat hail. Turbulent Mixing in Super- 
sonic Jets. PhD thesis, Penn State University, 

1991. 

8. C. C. Chang and C. Y. Kuo. Instability of a 
supersonic vortex sheet inside a circular duct. 
Phys. Fluids A, 5:2217 2228. 1993. 

9. K. Viswanathan, P. J. Morris, and G. Chen. 
Linear instability waves in supersonic jets con- 
fined in circular and non-circular ducts. J. 
Sound Vib., 171:231 253. 1994. 

10. C. K. W. Tam and D. E. Burton. Sound gener- 
ated by instability waves of supersonic flows, 
part 2. axisymmetric jets. J. Fluid Meek., 
138:273-295,* 1984. 

11. M. D. Dahl. The Aeroacoustics of Supersonic 
Coaxial Jets. PhD thesis, Penn State Univer- 
sity, 1994. 

12. A. Michalke and G. Hermann. On the inviscid 
instabilitv of a circular jet with external flow. 
J. Fluid Meek 114:343-359. 1982. 

13. M. E. Goldstein. Aeroacoustics. McGraw-Hill, 
Now York, 1976. 


NASA TM- 107441 


11 



REPORT DOCUMENTATION PAGE 


Form Approved 
OMB No. 0704-0188 


Public reporting burden for ttiis collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, 
gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect of this 
collection of information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for Information Operations and Reports, 1215 Jefferson 
Davis Highway, Suite 1204, Arlington, VA 22202-4302, and to the Office of Management and Budget, Paperwork Reduction Project (0704-0188), Washington, DC 20503. 


1 . AGENCY USE ONLY ( Leave blank) 


4. TITLE AND SUBTITLE 


2. REPORT DATE 

May 1997 


3. REPORT TYPE AND DATES COVERED 

Technical Memorandum 


5. FUNDING NUMBERS 


Effects of Acoustically Lined Cylindrical Ducts on Instability Waves in Confined 
Supersonic Jets 



WU-505-62-52 


7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135-3191 

8. PERFORMING ORGANIZATION 
REPORT NUMBER 

E- 10709 

9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 

10. SPONSORING/MONITORING 
AGENCY REPORT NUMBER 

National Aeronautics and Space Administration 
Washington, DC 20546-0001 

NASA TM- 107441 
AIAA-97-1600 

11. SUPPLEMENTARY NOTES 

Prepared for the Third Aeroacoustics Conference sponsored by the American Institute of Aeronautics and Astronautics and 
the Confederation of European Aerospace Societies, Atlanta, Georgia, May 12-14, 1997. Responsible person, Milo D. 
Dahl, organization code 2660, (216) 433-3578. 

12a. DISTRIBUTION/A VAILABILITY STATEMENT 

12b. DISTRIBUTION CODE 

Unclassified - Unlimited 
Subject Category 34 


This publication is available from the NASA Center for AeroSpace Information, (301) 621-0390. 



13. ABSTRACT (Maximum 200 words) 

The pressure disturbances generated by an instability wave in the shear layer of a supersonic jet are studied for an 
axi symmetric jet inside a lined cylindrical duct. For the supersonic jet, locally linear stability analysis with duct wall 
boundary conditions is used to calculate the eigenvalues and the eigenfunctions. These values are used to determine the 
growth rates and phase velocities of the instability waves and the radial pressure disturbance patterns. The study is 
confined to the dominant Kelvin-Helmholtz instability mode and to the region just downstream of the nozzle exit where 
the shear layer is growing but is still small in size compared to the radius of the duct. Numerical results are used to study 
the effects of changes in the outer flow, growth in the shear layer thickness, wall distance, wall impedance, and frequency. 
Results indicate that the effects of the duct wall on shear layer growth rates diminish as the outer flow increases. Also, 
wall reflections cause variations in growth rates depending on wall height and Strouhal number. These variations are due 
to the phase relationship between the outgoing and the reflected incoming pressure disturbances at the shear layer. The 
growth rate variations can be reduced and the maximum growth rate minimized by keeping the imaginary part of the 
impedance negative. 


14. SUBJECT TERMS 



Supersonic jets; Acoustic treatment; Instability waves; Ejectors; Confined jets 

17. SECURITY CLASSIFICATION 

18. SECURITY CLASSIFICATION 

19. SECURITY CLASSIFICATION 

OF REPORT 

OF THIS PAGE 

OF ABSTRACT 

Unclassified 

Unclassified 

Unclassified 


15. NUMBER OF PAGES 

13 


16. PRICE CODE 

A03 


20. LIMITATION OF ABSTRACT 


NSN 7540-01-280-5500 


Standard Form 298 (Rev. 2-89) 
Prescribed by ANSI Std. Z39-18 
298-102 























